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Abstract. Populations of relativistic electrons in clusters of galaxies should imprint a characteristic spectral dis- 
tortion signature on the Cosmic Microwave Background (CMB) photons - the relativistic Sunyaev-Zeldovich (SZ) 
effect. We investigate how sensitive PLANCK and successor experiments will be in detecting or constraining the 
total optical depth r rc i of relativistic electrons. We find an expected sensitivity of PLANCK of oy r<> , ~ 2 ■ 10 -5 
for an individual cluster, which is too insensitive to detect even optimistic scenarios, which predict r rc i ~ 10 -6 
to 10 -7 . However, by stacking the PLANCK-signal of all SZ detectable clusters the relativistic SZ effect can be 
statistically probed to a sensitivity level of oy l(s , ~ 10~ 7 . PLANCK successor experiments will be able to probe 
even conservative scenarios of relativistic electron populations in clusters of galaxies. 
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1. Introduction 
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Clusters of galaxies contain a several keV hot thermal 
plasma, which is able to Comptonize the CMB while 
its photons are traversing the cluster. The resulting ob- 
servable spectral distortions of the CMB - the SZ effect 
(Sunyaev & Zeldovich 1972) - are discriminated by the 
different nature of the energy forms of the scattering elec- 
trons: the thermal energy manifests itself in the thermal 
SZ (tSZ) effect, and the kinetic energy of a gas bulk motion 
in the kinematic SZ (kSZ) effect. The spectral characteris- 
tics of the tSZ and the kSZ effects are independent of the 
electrons' energies as long as the electrons' velocities are 
much smaller than the speed of light c. However, the ther- 
mal electrons in hot galaxy clusters with temperatures of 
k 10 keV have velocities of the order of 0.2 c, and there- 
fore relativistic corrections to the canonical tSZ spectral 
distortions have to be taken into account (Wright 1979). 

In addition to the mildly relativistic thermal electrons 
there are also highly relativistic electrons present in galaxy 
clusters, which have their own characteristic relativistic SZ 
(rSZ) spectral distortion signature (Enfilin & Kaiser 2000). 
It is the aim of this work to investigate how future sen- 
sitive spectral measurements of CMB distortions caused 
by galaxy clusters can be used in order to investigate or 
constrain such electron populations. 

We know about the existence of highly relativistic elec- 
trons in galaxy clusters since extended, cluster-wide radio 



emission, the so called cluster radio halos, has been de- 
tected in roughly one third of the X-ray brightest galaxy 
clusters (Giovannini et al. 1999). The radio halo emis- 
sion is produced by radio synchrotron emitting ultra- 
relativistic electrons with energies of the order of 10 GeV. 
Relativistic electrons at energies far below a GeV have es- 
caped our measurements so far since their synchrotron 
emission is at frequencies below the Earth ionosphere 
transparency cutoff, although their dim synchrotron self- 
Comptonized radiation should fall into our telescope's fre- 
quency bands (Enfilin & Sunyaev 2002). However, large 
numbers of such relativistic electrons are expected, not 
only as a low energy continuation of the observed radio 
emitting electron spectra, but also due to numerous phys- 
ical generation mechanisms. Such mechanisms are cluster 
merger shock waves, radio galaxy outflows, and electron 
production in inelastic interactions of relativistic protons 
with the thermal gas nucleons (Enfilin 1999b; Brunetti 
2002, for reviews). There are two very different habitats 
of relativistic electrons in galaxy clusters. The first is the 
relatively dense thermal gas, in which Coulomb cooling 
should quickly thermalize any significant electron popula- 
tion below 100 MeV (Sarazin 1999). The second are the 
fossil remnants of former radio galaxy cocoons, the radio 
quiet ghosts cavities (Enfilin 1999a), in which even MeV 
electrons may survive for cosmological times due to the 
low gas density there. The physical origin of any high en- 
ergy electrons is very interesting, and since our present 
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observational knowledge about these MeV to GeV elec- 
trons is poor (Enfilin & Biermann 1998, for a compila- 
tion of the available spectral information on Coma), any 
means to obtain additional information is highly desirable. 
In this paper we will analyze the ability of the PLANCK 
experiment or its successors to detect or constrain such 
electrons. 

In contrast to a thermal electron population 1 , the spec- 
tral shape of a relativistic electron population is poorly 
known. At GeV energies, they are well described by power 
laws. However their number density is generally domi- 
nated by the unknown lower end of their spectral distri- 
bution, which is expected to have a broad peak. Since the 
rSZ effect depends practically only on the total number of 
relativistic electrons and very little on their spectral distri- 
bution we approximate the spectrum by a delta function 
centered on that spectral peak. As will become clearer in 
the following, this idealized treatment is sufficient for our 
investigations of the sensitivity to detect the rSZ effect. 
A more realistic treatment would introduce several poorly 
constrained parameters like spectral index and low energy 
cutoff and its detailed shape, for which the natural choice 
is not obvious. 



2. Spectra of CMB distortions 

The PLANCK experiment should be able to identify 
~ 4 • 10 4 clusters of galaxies via their tSZ effect (e.g. 
Bartelmann 2001). In order to be able to do spectral pre- 
cision measurements of the CMB distortions caused by a 
galaxy cluster one has to suppress other contaminations as 
far as possible. A way to do this is optimal spatial filtering, 
which allows to remove efficiently foregrounds like galactic 
dust and synchrotron emission, but also CMB tempera- 
ture fluctuations, which have their spatial power mostly 
on large angular scales (Hachnclt & Tcgmark 1996). In the 
following we assume that such a filtering step was applied 
to the frequency maps of PLANCK, leading to galaxy clus- 
ter spectra which are mostly free of such contaminations. 
Since the galactic synchrotron and dust emission might 
have some power on small angular scales, which would 
partly survive such a filtering step, we further assume that 
the lowest and highest frequency channels of PLANCK are 
used to remove these residual contaminations. 
The CMB-blackbody spectrum 



I v = i i(x) 
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(1) 



1 The reported high energy X-ray excess of the Coma galaxy 
cluster (Fusco-Femiano et al. 1999) triggered some speculations 
about the existence of a supra-thermal electron pupulation in 
between 10 and 100 keV, which also would produce a unique 
SZ signature (Enfilin et al. 1999; Enfilin & Kaiser 2000; Blasi 
et al. 2000; Blasi 2000; Liang et al. 2002; Colafrancesco et al. 
2003). However, such a population is questioned on theoreti- 
cal grounds (Petrosian 2001), and even the high energy X-ray 
excess itself is under debate (Rossetti & Molendi 2004; Fusco- 
Femiano et al. 2004). 



with x = hu/kT and i = 2 (fcT cm b) 3 /(/ic) 2 becomes dis- 
torted due to the inverse Compton scattering off the elec- 
trons in a galaxy cluster. The spectral distortions are de- 
scribed by 

6i(x) = g(x) y gas (1 + S(x, T e )) - h(x) f3 gas r gas 

+(j(x,p)-i(x))T lel . (2) 

The first rhs term describes the tSZ distortions, which 
have a shape given by 

e x + 1 
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(e x - l) 2 V e x - 1 
and a magnitude described by the the Comptonization 
parameter y gas = a^/(m c c 2 ) Jdln cgas kT e . For non- 
relativistic electrons the relativistic correction term is 
zero, S(x,T e ) — 0, but for hot clusters even the thermal 
electrons are slightly relativistic, which will modify the 
thermal SZ effect (Wright 1979). These corrections are 
easily calculated (see e.g. Rephaeli 1995; Itoh & Nozawa 
2003; Enfilin & Kaiser 2000; Dolgov et al. 2001), and can 
be used to measure the cluster temperature purely from 
SZ observations (Hansen et al. 2002). For the PLANCK 
channels this is most easily implemented using the simple 
fitting formulae presented in Diego et al. (2003). 

The second rhs term in Eq. 2 gives the kSZ distortions, 
which have the spectral shape 



h(x) 
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and depend on /3 gas , the average line-of-sight streaming 
velocity of the thermal gas (v gas = /3gas c, (3 gas > if gas 
is approaching the observer), and the Thomson optical 
depth T gas = cr T Jdln Ctgas . 

Finally, the last rhs term in Eq. 2 describes the dis- 
tortion by relativistic electrons, which we approximate 
to be mono-energetically distributed. We characterize the 
electron energy via the dimensionlcss momentum p = 
P e / (m c c), where P e is the usual electron momentum. The 
optical depth r rc i = ctt Jdln e ^ e i determines the fraction 
of strongly scattered photons. Here, n c ^ IC \ is the number 
density of relativistic electrons. The spectral distortions 
Sirci = (j(x,p) — i(x))r ro i are described by two terms: 
The IC scattering redistributes the CMB photons, which 
means that photons disappear from the CMB spectrum 
according to — r IC \i(x), but reappear at different frequen- 
cies like T rc \j(p,x), which depends on the electron mo- 
mentum. In general, <5i ro i has to be calculated numeri- 
cally (Enfilin & Kaiser 2000), however the non- and ultra- 
relativistic limits can be treated analytically. 

In the non- relativistic limit (p <C 1), a mono-energetic 
electron population leads to spectral distortions which are 
identical to those of the well-known tSZ effect with 

„2 



<5«rcl = ^-g{x)T Ie l. 



(5) 



In the ultra-relativistic limit (p 3> 1) CMB photons simply 
disappear from the spectral range of interest since they are 
scattered to very high frequencies, leading to 

Si Ic i = -i(x) T re l. (6) 
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As long as all electrons are in one of these two limit- 
ing cases, the relative shape of the spectral distortion is 
completely independent of the details of the electron spec- 
trum, only the distortion strength is either proportional to 
the electrons total energy content (non-relativistic limit) 
or their number density (ultra-relativistic limit). 

Thus, only in the trans- relativistic regime (p ~ 1) 
the shape of the spectral distortions depends on de- 
tails of the electron spectrum. Since we are interested 
in the ultra-relativistic case, our approximation of mono- 
energetic electron spectra is fine. By keeping the full rSZ 
distortion term Si rc i = {j(x,p) — i(x))r rc i (for which we 
use the formulae given in Enfilin & Kaiser 2000), we make 
sure that our formalism tells us when we enter the trans- 
relativistic regime. For small p the rSZ spectral distortion 
approaches the one of the tSZ effect and thereby our sen- 
sitivity to discriminate relativistic electrons decreases. 

An order of magnitude estimate of the expected rSZ 
effect is in order. Any scattering of a CMB photon with 
an ultra-relativistic electron removes the photon from the 
typical CMB energy range. The relativistic electron op- 
tical depth T rc i can be expressed in terms of the tSZ y- 
parameter, if a characteristic ratio of the relativistic elec- 
tron to thermal gas energy density A^ei^ can be assumed 
to hold approximately on macroscopic scales throughout 
the cluster: 

r ro , = — — —, (7) 

where j c = E c /(m c c 2 ) = y/l + p 2 . Inserting some opti- 
mistic values for a Coma-like galaxy cluster like y gas ~ 
10~ 4 , X Te \ iB ~ 0.03, and 7 C <~ 10 for a relativistic elec- 
tron population residing in radio quiet ghost cavities 
(X Te i je ~ 0.1, and 7 C ~ 300 for relativistic electrons mixed 
with the thermal gas), one gets 5I le \/I = — r ro i ~ — 10~ 6 
(or 5I Te \/I ~ — 10 -7 respectively). 

3. Extracting the rSZ effect 

Our goal is to extract the rSZ spectral signature from 
PLANCK experiment multi-frequency measurements. In 
order to estimate the significance of our measurement, 
we apply a x 2 -statistics and search for the x 2 -minimum 
with respect to the unknown parameters. Minimization 
of the parameters within the ranges given in table 1 is 
done through an extention of the publicly available code 
sasz (Hansen 2003), based on the technique of simulated 
annealing. In our case the vector in parameter space, x, is 
a 5 dimensional vector, x = (y, T e , Vg asi T rc i, p). We choose 
a simple simulated annealing exponential cooling scheme, 
Tj = cTj-i, where c w 0.8 — 1.0, and we use T = 1 and 

Tfinal - 10- 12 . 



In fig. 1 we plot A% 2 as a function of the unknown r rc i 
for the PLANCK satellite, where we optimistically assume 
that only the first and last channels (30 and 857 GHz) are 



Parameter Allowed range 
y 10~' to 10"^ 
T e to 25 keV 

i; gas -5000 to 5000 km/sec 
r ro i -10~ 4 to 10~ 4 
p 0.01 to 1000 

Table 1. The allowed parameter range. 




T rel 



Fig. 1. A\ 2 as a function of r rc i, marginalized over the other 
parameters, (y, T e , w gas , p). la error-bar of r rc i is about 2-10~ 5 . 

needed for removal of cluster-scale foreground contamina- 
tion (like point sources, and galactic synchrotron and dust 
emission). We assume that the 'true' signal is T rcl = 10~ 6 
(the other unknown parameters are y = 10~ 5 , T e = 5 keV, 
Vgas = 600 km/sec, p — 200), and we maximize over all 
the unknown parameters. We find ler error-bars for the 
relativistic optical depth r re i just about 0> rol « 2 • 10~ 5 . If 
one instead excludes also the 545 GHz channel, then the 
error-bars are increased by approximately a factor of 2. It 
is comforting to note that the best fit point is indeed very 
near the 'truth' 10~ 6 . 

We next have to test for parameter degeneracy. As is 
well known, parameter degeneracy is very important for 
the error-bars on the normal SZ parameters. E.g. a large 
negative peculiar velocity can mimic a higher tempera- 
ture (and slightly changed Compton parameter). In this 
way a large negative peculiar velocity of a cluster results 
in larger error-bars on the determination of the peculiar 
velocity than in the case of a large positive peculiar ve- 
locity (Aghanim et al. 2003). We ran a sample of different 
parameters (T e = 2 - 10 keV, w gas = -1000 - +1000 
km/sec, y = 10~ 6 — 10~ 4 ), and we find both very little 
parameter degeneracy with r re i, and that the error-bar on 
r ro i remains 2 • 10~ 5 within a factor of 2 for all the clusters 
in this sample. This also implies that systematic shifts of 
the other cluster paramters due to non-removable contam- 
ination should not strongly alter our conclusions. 

As can be clearly seen on figure 2 the determination 
of T re i gives much larger error-bars for p < 3 (as explained 
in section 2), compared to the case p > 3 with small 
and constant error-bars. t ic \ is therefore practically non- 
determinable forp < 1. However, the unproblcmatic p > 3 
case is exactly the range we are interested in. 

We therefore conclude that the error-bar on the deter- 
mination of Trci will be about 2- 10~ 5 for any single typical 
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Fig. 2. The error-bar on r ro i as a function of the dimension less 
electron momentum, p. The error-bar on T ro i is independent of 
p as long as p > 3. 

clusters detectable by PLANCK, roughly independent of 
other cluster parameters. 

4. Conclusion 

We estimate the sensitivity of PLANCK to detect the op- 
tical depth r re i of relativistic electrons to be ay rel ~ 2- 10~ 5 
for a single cluster. This is too insensitive to detect even 
optimistic senarios for the relativistic electron population, 
which give values of r rc i ~ 1CT 7 to 10~ 6 . 

However, there is no significant parameter degeneracy 
with other cluster parameters for the measurement of the 
optical depth t ic \ of relativistic electrons. This allows one 
to stack the weak signals of the large number of galaxy 
clusters to be detected by PLANCK with the thermal SZ 
effect. Since PLANCK is expected to observe of the order 
of N c \ w 4 • 10 4 individual clusters (e.g. Bartelmann & 
White 2002), the expected sensitivity for the statistical 
detection of the relativistic SZ signal should be of the order 
oi tat = a T , N^ 1 ^ ~ 10~ 7 . This is sufficient to confirm or 

' rcl rcl CI 

refute optimistic scenarios. 

Future experiments with increased sensitivities may be 
able to test for even smaller optical depth of relativistic 
electron populations. An experiment with 4 observing fre- 
quencies at 90, 150, 217 and 270 GHz, and with South- 
Pole- Telescope 2 -like sensitivity of 0.1 [iK would allow to 
probe relativsitic electrons down to <7 Trel ~ 2 • 10~ 6 for 
an individual cluster. Thus optimistic theoretical scenar- 
ios could be tested using a few individual clusters, while 
statistical measurements would probe even deep into con- 
servative scenarios of relativistic electron populations in 
galaxy clusters. 

We conclude that sensitive CMB experiments like 
PLANCK allow to constrain or examine the presently 
poorly known low energy part of the existing relativis- 
tic electron populations in galaxy clusters. Therefore they 
provide a valuable vehicle to explore this terra incognita 
in electron energies in between MeV and GeV. 
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